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• We present a new class of 3D black hole initial data sets for numerical relativity. These 

' data sets go beyond the axisymmetric, "gravity wave plus rotating black hole" single 

black hole data sets by creating a dynamic, distorted hole with adjustable distortion pa- 
rameters in 3D. These data sets extend our existing test beds for 3D numerical relativity, 
, representing the late stages of binary black hole collisions resulting from on-axis collision 

or 3D spiralling coalescence, and should provide insight into the physics of such systems. 
We describe the construction of these sets, the properties for a number of example cases, 
and report on progress evolving them. 

1 Introduction 

\^ ■ In all that follows, the standard ADM 3+1 decomposition is usedB along with the 

standard York decomposition of the initial SfJue problcma This work extends 
previous studies of axisymmetric black holescru to allow for a dependence on the 
azimuthal coordinate. The basic idea is to combine the pure gravitational wave 
construction due to Brill and combine it with a black hole to form adjustable, 
| highly distorted black hole initial data sets. 

U" 1 2 Initial Value Problem I: Non- Axisymmetric Distorted Black Holes 

The first type of solution we consider for a non-axisymmetric gravity wave is the 
Brill wave, form superimposed on an Einstein- Rosen bridge, used by 2D simulations 
at NCSAtl In these time-symmetric spacetimes, the three metric is given by 



X 



^ ; ds 2 = i/j 4 (e 2q [dr) 2 + dO 2 } + sin 2 6 d(/> 2 ) . (1) 

In this equation, q may be chosen freely (so long as the appropriate boundary 
conditions are respected), but the form we use here is similar to that used by 
previous NCSA worm 

q = Q o sin" 9 (l + c cos 2 <j>)g (2) 
g = e -(v-vo) 2 /<? 2 + e ~(ri+no) 2 /<r 2 _ (3) 

In the above n is an even integer, constrained such that n > 2, which specifies the 
angular dependence of the Brill wave. The parameter a specifies its width, and 770 
its radial position. The parameter c specifies the non-axisymmetric aspect of the 
wave, and when c = we recover the axisymmetric caseo. Once these parameters 
are supplied, we numerically solve the Hamiltonian constraint equation for ■0, the 
conformal factor. The initial value code was tested against the results of the 2D 
code for accuracy, and convergence studies show that it converges to second order. 
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3 Initial Value Problem II: Non-Axisymmetric Distorted Black Holes 
with Angular Momentum 



Another form of distorted black hole initial data, one which allows us to have both a 
non-axisymmctric gravity wave and angular momentum, is based on a conformally 
flat space (q — 0). Here, the gravity wave is constructed by specifying the conformal 
extrinsic curvature. Its parameters are an angular momentum J, a (^-dependence 
given by m (normally we use either m — 4 so that the IVP is isometric under the 
operation </> — > </>±7r/2, or we use m = to facilitate comparison with 2D codes), and 
a function of r\ that is typically chosen to be the same as g above, and is arbitrary 
so long as it vanishes as rj — > oo and is symmetric at 77 = 0. Here we sketch the 



procedure of constructing the solution to the constraints: writing 
take 



Ki. 





h 


hi 















m sin(m0) 


hi 


h x — ho 





+ 


















—h x sin 2 9 




. K n $ 








%l> Kij we 



(4) 



K^ — sin 9 ( 3 J + cos (m 



1 



-da (sin 



Kg 



= cos (m(f>) sin 9 (m 2 v + sin 2 9d v g) 
If we now make the substitutions: 
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then the 77- and 0-momentum constraints become: 

{dl-l + d 2 B +2cotedg)n = - (4 cos 9 + sin 9d e ) d^g 



(d 2 -l + d e 2 + 2cotedg)A 



cscedd,v — dgh x 



sin edrjdrj, 



(5) 
(6) 

(7) 
(8) 



(9) 
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Carefully choosing the forms and boundary values of the functions on the right hand 
side will then allow us to solve the constraints analytically in terms of the function 
g and its derivatives. Details will be published in a future paper. 



4 Evolutions of type I initial data 

To evolve this data, we first interpolate it from the spheroidal (77, 6, (j)) coordinate 
system to a Cartesian coordinate system used by our evolution codes. For the 
present we limit our evolutions to a grid with an outer boundary placed at 18M, 
so the results we obtain should not be taken as more than a feasibility study. 
Higher resolution runs, on much larger grids, are in progress and will be publishedj'Q 
elsewhere. Nevertheless, we are able to evolve both axisymmetric data (c = 0), and 
non-axisymmctric data (c ^ 0) and obtain qualitatively similar results. The metric 
elements evolve in the same fashion, with the same qualitative features we are 



2 



familiar with from 2D codes. Below we show two figures. In Figure I we see the 
growth of the metric function for a run parameterized by (Qo = —0.1, 770 = 0.0, a = 
1.0, n = 4, c = 0.5). The growth of the metric function here is familiar from all ID 
and 2D studies performed by the group. In Figure II we see the beginnings of a 
waveform extracted at several radii and can watch its progress as it travels outward. 
The wavelength matches well with the value predicted from linear theory, despite 
the nearness of the outer boundary. 
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Figure 1: Metric function g xx evaluated 
along the x-axis and plotted at several 
times. 




Time (M ADM ) 

Figure 2: The I = 2 Zerilli function ex- 
tracted at three different radii and plot- 
ted as a function of time. The wave 
function is clearly propagating outward. 



5 Future Work 

As of this writing, much of this work has already been extended!. Furthermore, 
in the near future we plan to publish full descriptions of these data sets, their 
implementations, and to describe in detail comparisons with 2D code results (where 
the outer boundary may be placed arbitrarily far away), and convergence tests. 
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